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■ These lectures give an overview of aspects of missing Ej- signatures from new 
' physics at the LHC, along with their important connection to dark matter 
i physics. Mostly, I will concentrate on supersymmetric (SUSY) sources of Ej>, 

but will also mention Little Higgs models with T-parity (LHT) and univer- 
sal extra dimensions (UED) models with K iif-parity. Lecture 1 covers SUSY 
basics, model building and spectra computation. Lecture 2 addresses sparticle 
production and decay mechanisms at hadron colliders and event generation. 
Lecture 3 covers SUSY signatures at LHC, along with LHT and UED signa- 
tures for comparison. In Lecture 4, I address the dark matter connection, and 
how direct and indirect dark matter searches, along with LHC collider searches, 
may allow us to both discover and characterize dark matter in the next several 
years. Finally, the interesting scenario of Yukawa-unified SUSY is examined; 

■ this case works best if the dark matter turns out to be a mixture of axion/axino 
' states, rather than neutralinos. 
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1. Introduction: Et collider signatures and the dark matter 
connection 

I have been assigned the topic of discussing missing transverse energy sig- 
natures at hadron colliders, especially the CERN LHC, which should be- 
gin producing data in 2009. Every collider event- be they lepton-lepton, 
lepton-hadron or hadron-hadron interactions- will contain some amount of 
missing energy just due to the fact that energy measurements are not 
perfectly precise. Hadronic events will usually contain more Ifi than leptonic 
events since the energy resolution of hadronic calorimeters is not as precise 
as is electromagnetic resolution. In addition, collider detectors are imper- 
fect devices, often containing un-instrumented regions, cracks, hot or dead 
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calorimeter cells, and only partial coverage of the 4tt steradians surround- 
ing the interaction region: always, an allowance must be made at least to 
allow the beam pipe to enter the detector. Also, missing energy can arise 
from cosmic rays entering the detector, or beam-gas collisions, multiple 
scattering etc. 

At e + e~ colliders, $ can be directly measured since the incoming beam 
energies are well-known (aside from bremsstrahlung/beamstrahlung effects 
at very high energy e + e~ colliders) . At pp or pp colliders, the hard scattering 
events which are likely to produce new physics are initiated by quark-quark, 
quark-gluon, or gluon-gluon collisions, and the partonic constituents of the 
proton carry only a fraction of the beam energy. The proton constituents not 
participating in the hard scattering will carry an unknown fraction of beam 
energy into the un-instrumented or poorly instrumented forward region, so 
that for hadron colliders, only missing transverse energy ($t) is meaningful. 
In these lectures, we will concentrate on pp colliders, in anticipation of the 
first forthcoming collisions at the CERN LHC. 

Along with these imperfect detector effects, fir can arise in Standard 
Model (SM) production processes wherein neutrinos are produced, either 
directly or via particle decays. For instance, the key signature for W boson 
production at colliders was the presence of a hard electron or muon balanced 
by missing energy coming from escaping neutrinos produced in W~ — > l~vi 
(or charge conjugate reaction) decay (here, I — e, \i or r). In addition, 
neutrinos will be produced by semi-leptonic decays of heavy flavors and rs, 
and indeed $t was an integral component of tt production events followed 
by t — > bW decays at the Fermilab Tevatron, and led to discovery of the 
top quark. 

In these lectures, we will focus mainly on new physics reactions which 
lead to events containing large amounts of $r- The major motivation nowa- 
days that $t signatures from new physics should appear at LHC comes from 
cosmology. A large array of data- coming from measurements of galactic 
rotation curves, galaxy cluster velocity profiles, gravitational lensing, hot 
gas in galaxy clusters, light element abundances in light of Big Bang nu- 
cleosynthesis, large scale structure simulations and measurements and es- 
pecially measurements of anisotropies in the cosmic microwave background 
radiation (CMB)- all point to a consistent picture of a universe which is 
constructed of 

• baryons: ~ 4% 

• dark matter: <~ 21% 

• dark energy: ~ 75% 
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• vs, 7s: a tiny fraction. 

Thus, particles present in the SM of particle physics constitute only <~ 4% 
of the universe's energy budget, while the main portion is comprised of 
dark energy (DE) which causes an accelerating expansion of the universe 
and dark matter (DM). There is an ongoing program of experiments under 
way or under development which will probe the dark energy. The goal is 
to try to tell if DE comes from Einstein's cosmological constant (which 
theoretical prejudice says really ought to be there), or something quite 
different. Measurements focus on the parameters entering the dark energy 
equation of state. 

For dark matter, there exist good reasons to believe that it will be 
connected to new physics arising at the weak scale: exactly the energy 
regime to be probed by the LHC. The density of dark matter is becoming 
known with increasing precision. The latest WMAP5 measurements 1 find 

VcDMh 2 = 0.110 ±0.006 (1) 

where Sic dm is the cold DM density divided by the critical closure density 
of the universe p c , and h is the scaled Hubble constant: h = 0.73^'^- While 
the cosmic DM density is well-known, the identity of the DM particle(s) is 
completely unknown. Nevertheless, we do know some properties of the dark 
matter: it must have mass, it must be electrically neutral (non-interacting 
with light) and likely color neutral, and it must have a non-relativistic 
velocity profile, i.e. be what is termed cold dark matter (CDM). The reason 
it must be cold is that it must seed structure formation in the universe, i.e. 
it must clump. That means its velocity must be able to drop below its escape 
velocity; otherwise, it would disperse, and not seed structure formation. We 
already know of one form of dark matter: cosmic neutrinos. However, the 
SM neutrinos are so light that they must constitute hot dark matter. 

The theoretical literature contains many possible candidate DM parti- 
cles. A plot of some of these is shown in Fig. 1 (adapted from L. Roszkowski) 
in the mass vs. interaction cross section plane. Some of these candidates 
emerge from attempts to solve longstanding problems in particle physics. 
For instance, the axion arises from the Peccei-Quinn solution to the strong 
CP problem. The weakly interacting massive particles (WIMPs) frequently 
arise from models which attempt to explain the origin of electroweak sym- 
metry breaking. WIMP dark matter candidates include the lightest neu- 
tralino of models with weak scale supersymmetry, while Kaluza-Klein pho- 
tons arise in models with universal extra dimensions, and lightest T-odd 
particles arise in Little Higgs models with a conserved T-parity. 
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Fig. 1. Dark matter candidates in the mass versus interaction strength plane, 
taken from Ref.: z http://www.science.doe.gov/hep/hepap_reports.shtm. See also, 
L. Roszkowski, Ref. 3 



WIMP particles have an additional motivation, which has recently been 
coined the "WIMP miracle". 4 The idea here is to assume the existence 
of a dark matter particle which was once in thermal equilibrium at high 
temperatures in the early universe. The equilibrium abundance is easily 
calculated from thermodynamics and provides one boundary condition. As 
the universe expands and cools, the DM particle will drop out of thermal 
cquilbrium (freeze-out); while DM will no longer be produced thermally, 
it can still annihilate away. The subsequent time evolution of the dark 
matter density is governed by the Boltzmann equation as formulated for 
the Friedman-Robertson- Walker universe: 



dn/dt = —3Hn ~ (crv re i)(n 2 



(2) 



wherein the —3Hn terms gives rise to dilution due to the Hubble expansion, 
and the (av re i) term involves the dark matter annihilation cross section. 
The Boltzmann equation is solved in various textbooks, 5 and the solution 
is given approximately by 



ilh 2 = 
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p c /h 2 \ng* ) Mpi (av) ' 

where Sq is the current entropy density, is the effective degrees of freedom 
and Xf is the freeze-out temperature ~ The quantities entering ths 
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expression are unrelated to one another and several have large exponents. 
Filling in numerical values, they conspire to give 

^ ~ °^J* (4) 

(av) 

so that in order to obtain the measured dark matter abundance flh 2 ~ 0.11, 
we expect (av) to be a weak scale sized cross section).. Inputing a typical 
weak scale cross section, 

a= W (5) 
and setting a ~ 1 pb, we find m ~ 100 GeV! The lesson here is an important 
one: that the measured dark matter density is consistent with the annihi- 
lation cross section expected of a weakly interacting particle with mass of 
order the weak scale. This auspicious relation amongst dark matter density, 
annihilation cross section and particle mass is sometimes refered to as the 
WIMP miracle, and seemingly offers independent astrophysical evidence for 
new physics at or around the weak scale: in this case the new physics is 
associated with the cosmic density of dark matter. If such a WIMP particle 
exists, then it should be produced at large rates at the CERN LHC. Since 
it is dark matter, once produced in a collider experiment, it would escape 
normal detection. If produced in association with other SM particles (such 
as quarks and gluons), then it will carry away missing (transverse) energy, 
and will give rise to new physics signatures including jets plus JSt- It is 
this connection of weak scale physics and dark matter which makes the jets 
+ |?t signature so important and exciting at the LHC. 

2. Lecture 1: SUSY basics: model building and spectra 
calculation 

2.1. The Wess-Zumino model 

We begin with a toy model due to Wess and Zumino, from 1974. 6 The 
model Lagrangian is given by 

£ ^kin. "1~ ^mass (6) 

with 

Adn. = \{3^Af + \{d mu Bf + \{F 2 + G 2 ) 



£ mass = -m[^V -GA- FB], 
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where A and B are real scalar fields, ip is a 4-component Majorana spinor 
field and F and G are non-propagating auxiliary fields. The auxiliary fields 
may be eliminated from the Lagrangian via the Euler-Lagrange equations: 
F = —mB and G = —mA. 

But instead, impose the transformations 

A—>A + 5A with 5 A = iaj 5 ip, 
SB = —clip, 

Sip = —Fa + iG^a+ cf^Act + i tfBa, 
SF = ia (pip, 
SG = 0:75 cf^). 

Note that the transformation mixes bosonic and fermionic fields. It is called 
a super symmetry transformation. Using Majorana bilinear re-arrangements 
(such as V>x = — xVO) the product rule • g) = d^f -g + f- d^g and some 
algebra, we find that C — > C + SC with 

1 i i 1 

<5£ kin = ^{--a^n $Bip + -0:757^ $M + 2 F ^^ + T^" 757 ^)' 

^ mass = ^(772^0:75 7^ V + imBa-f^ip). 

Since the Lagrangian changes by a total derivative, the action S = J Cd 4 x is 
invariant! (owing to Gauss' law in 4-d j v d^xd^A^ = J gv daA^n^). Under 
supersymmetry transformations, the Lagrangian always transforms into a 
total derivative, leaving the action invariant. 

The Wess-Zumino (WZ) model as presented in 1974 exhibited a quali- 
tatively new type of symmetry transformation with astonishing properties: 
one being that the usual quantum field theory quadratic divergences associ- 
ated with scalar fields all cancel! Theories with this sort of property might 
allow for a solution of the so-called "gauge hierarchy problem" of grand 
unified theories (GUTs). 

2.2. Superfield formalism 

While the WZ model was very compelling, it was constructed by hand, 
and offered no insight into how to construct general supcrsymmetric mod- 
els. A short while later, Salam and Strathdee 7 developed the superfield 
formalism, which did provide rules for general SUSY model building. 
They first expanded the usual four-component spacetime x M to super space: 
x^ — ► (x M , 9 a ), where 9 a (a = 1 — 4) represents four anti- commuting (Grass- 
mann vaued) dimensions arranged as a Majorana spinor. Then they intro- 
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duced general superfields: 

0)=S- - l -{hsO)M + \{96)N + l(h>1vO)V 

left chiral scalar superfields (LCSSF): 

S L (x, 6) = S(x) + iV20ip L (x) + m L T{x) (7) 
(where x ^ = x^ + 5^757^^) and right chiral scalar superfields (RCSSF): 

S R (x, 9) = 5(xt) - iV26ip R (tf) - m R T{x^). (8) 
The multiplication table for superfields is as follows: 

• LCSSF x LCSSF = LCSSF 

• RCSSF x RCSSF = RCSSF 

• LCSSF x RCSSF = a general superfield. 

Moreover, under a supersymmetry transformation, the D-term of a general 
superfield and the F term of a LCSSF or a RCSSF always transforms into 
a total derivative. 

This last observation gave the key insight allowing for construction of 
general supersymmetric models. The superpotential / is a function of LC- 
SSFs only, and hence by the multiplication table above is itself a LCSSF. 
Its F term transforms into a total derivative under SUSY, and hence is a 
candidate Lagrangian. Also, the Kahler potential K is introduced: it is a 
function of products of RCSSFs times LCSSFs: it is therefore a general su- 
perfield and its D-term transforms into a total derivative under SUSY, and 
hence its D-term is a candidate Lagrangian. The choice K = fields 
is completely general for a renormalizable theory. Augmenting the above 
superfields with gauge superfields (ones containing spin-1 gauge bosons and 
spin-i gauginos all transforming under the adjoint representation) allows 
one to construct locally gauge invariant, renormalizable field theories which 
are also invariant under SUSY transformations. 

A Master Lagrangian can then be constructed. Its form is given in Ch. 
5 of Ref. 8 It allows for a recipe for SUSY model building: 

• select the gauge group (simple or direct product), 

• select the matter and Higgs representations, 

• write down a superpotential (gauge invariant and renormalizable), 

• plug all terms into the Master Lagrangian, calculate the physical 
states and Feynman rules, and you are good to go! 



January 29, 2009 20:12 



WSPC - Proceedings Trim Size: 9in x 6in tasi 



8 

2.3. Minimal Super symmetric Standard Model (MSSM) 

The superficld formalism 8-10 allows for the construction of a supcrsymmet- 
ric version of the Standard Model, known as the Minimal Supersymmetric 
Standard Model, or MSSM. For each quark and lepton of the SM, the 
MSSM necessarily includes spin-0 superpartners q~L and qr along with £l 
and £r, whose gauge quantum numbers are fixed to be the known gauge 
quantum numbers of the corresponding SM fcrmions. Thus, for example, 
the right-handed up quark scalar (usually denoted by ur) is a color-triplet, 
weak isosinglet with the same weak hypercharge 4/3 as the right-handed 
up-quark. The MSSM thus includes a plethora of new scalar states: &l, 
Cfl, v e L, ul, ur, dj,, da in the first generation, together with analogous 
states for the other two generations. Spin-zero squark partners of quarks 
with large Yukawa couplings undergo left-right mixing: thus, the tj, and tn 
states mix to form mass eigenstates - t\ and I2 - ordered from lowest to 
highest mass. 

The spin-0 Higgs bosons are embedded in Higgs supcrficlds, so that the 
MSSM also includes spin-^ higgsinos. Unlike in the SM, the same Higgs 
doublet cannot give a mass to both up- and down- type fermions without 
catastrophically breaking the underlying supersymmetry. Thus the MSSM 
includes two Higgs doublets instead of one as in the SM. This gives rise to 
a richer spectrum of physical Higgs particles, including neutral light h and 
heavy H scalars, a pseudoscalar A and a pair of charged Higgs bosons H ± . 

The gauge sector of the MSSM contains gauge bosons along with spin- 
half gauginos in the adjoint representation of the gauge group: thus, along 
with eight colored gluons, the MSSM contains eight colored spin-i gluinos. 
Upon electroweak symmetry breaking, the four gauginos of SU (2)l x U(1)y 
mix (just as the SU(2)l and U(1)y gauge bosons mix) amongst themselves 
and the higgsinos, to form charginos - and - and neutralinos - Z\ , 
Z 2 , Zy, and Z^. The Z\ state, the lightest neutralino, is often the lightest 
supersymmetric particle (LSP), and turns out to be an excellent WIMP 
candidate for CDM in the universe. 

If nature is perfectly supersymmetric, then the spin-0 superpartners 
would have exactly the same mass as the corresponding fermions. Charged 
spin-0 partners of the electron with a mass of 0.51 MeV could not have 
evaded experimental detection. Their non-observation leads us to conclude 
that SUSY must be a broken symmetry. In the MSSM, SUSY is broken 
explicitly by including so-called soft SUSY breaking (SSB) terms in the 
Lagrangian. The SSB terms preserve the desirable features of SUSY, such 
as the stabilization of the scalar sector in the presence of radiative correc- 
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tions, while lifting the superpartner masses in accord with what is necessary 
from experiment. It is important to note that the equality of dimension- 
less couplings between particles and their superpartners is still preserved 
(modulo small effects of radiative corrections): in particular, phenomenolog- 
ically important gauge interactions of superpartners and the corresponding 
interactions of gauginos remain (largely) unaffected by the SSB terms. 

2.4. Supergravity 

The addition of the SSB Lagrangian terms may seem ad-hoc and ugly. It 
would be elegant if instead supcrsymmetry could be spontaneously broken. 
But it was recognized in the early to mid- 1 980 's that models where global 
SUSY is spontaneously broken at the weak scale ran into serious difficulties. 
The situation is very different if we elevate SUSY from a global symmetry 
to a local one. In local SUSY, we are forced to include the graviton/gravitino 
super-multiplet into the theory, in much the same way that we have to in- 
clude spin-1 gauge fields to maintain local gauge invariance of Yang-Mills 
theories. Theories with local SUSY are known as supergravity (SUGRA) 
theories because they are supersymmetric and necessarily include gravity. 
Moreover, the gravitational sector of the theory reduces to general relativ- 
ity in the classical limit. Within the framework of SUGRA, it is possible 
to add an additional sector whose dynamics spontaneously breaks SUSY 
but which interacts with SM particles and their superpartners only via 
gravity (the so-called hidden sector). The spontaneous breakdown of su- 
persymmetry results in a mass for the gravitino in the same way that in 
local gauge theories gauge bosons acquire mass by the Higgs mechanism. 
This is, therefore, referred to as the super-Higgs mechanism. The remark- 
able thing is that because of the gravitational couplng between the hidden 
and the MSSM sectors, the effects of spontaneous supersymmetry breaking 
in the hidden sector are conveyed to the MSSM sector, and (provided the 
SUSY-breaking scale in the hidden sector is appropriately chosen) weak 
scale SSB terms that lift the undesirable degeneracies between the masses 
of SM particles and their superpartners are automatically induced. Indeed, 
in the limit where Mpi — > oo (keeping the gravitino mass fixed), we recover 
a global SUSY theory along with the desired SSB terms! The gravitino 
typically has a weak scale mass and in many cases decouples from collider 
experiments because of its tiny gravitational couplings. For reasons that 
we cannot discuss here, these locally supersymmetric models are free 8 " 10 of 
many of the difficulties that plague globally supersymmetric models. 
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2.5. SUGRA GUTs 

Motivated by the successful unification of gauge couplings at a scale 
Mqut ~ 2 x f0 16 GcV in the MSSM, we are led to construct a grand 
unified theory (GUT) based on local supersymmetry. In this case, the the- 
ory renormalized at Q = Mgu T contains just one gaugino mass parameter 
mi/2- Renormalization effects then split the physical gaugino masses in the 
same way the measured values of the gauge couplings arise from a single 
unified GUT scale gauge coupling. In general, supergravity models give rise 
to complicated mass matrices for the scalar superpartners of quarks and 
leptons, with concomitant flavor violation beyond acceptable levels. How- 
ever, in models with universal soft SUSY breaking terms, a super-GIM 
mechanism suppresses flavor violating processes. 11 In what has come to be 
known as the minimal supergravity (mSUGRA) model, a universal scalar 
mass mo and also a universal SSB scalar coupling Aq are assumed to exist 
at a high scale Q <~ Mqut — Mp\. The physical masses of squarks and 
slcptons are split after renormalization, and can be calculated using renor- 
malization group techniques. Typically, in the mSUGRA model, we have 
m,q ~ rrij; L ~ m i R - Although the Higgs scalar mass parameters also start 
off at the common value m at the high scale, the large value of the top 
quark Yukawa coupling drives the corresponding squared mass parameter 
to negative values and EWSB is radiatively broken. Within this framework, 
the masses and couplings required for phenomenology are fixed by just a 
handful of parameters which are usually taken to be, 

to , m 1 / 2 , A a , tan/5, and sign(fi). (9) 

Here, tan /3 is the ratio of the vacuum expectation values of the Higgs fields 
that give masses to up and down type fermions, and fj, is the supersym- 
metric higgsino mass parameter whose magnitude is fixed to reproduce the 
measured value of Mz- If all parameters are real, then potentially large 
CP-violating effects are suppressed as well. Computer codes such as Isajet, 
SuSpcct, SoftSUSY and Spheno that calculate the full spectrum of sparticle 
and Higgs boson masses are publicly available. 12 

2.6. Some realistic SUSY models 

The mSUGRA model (sometimes referred to as the constrained MSSM or 
CMSSM) serves as a paradigm for many SUSY phenomenological analy- 
ses. However, it is important to remember that it is based on many as- 
sumptions that can be tested in future collider experiments but which 
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may prove to be incorrect. For instance, in many GUT theories, it is com- 
mon to get non-universal SSB parameters. In addition, there are other 
messenger mechanisms besides gravity. In gauge-mediated SUSY break- 
ing models (GMSB), 13 a special messenger sector is included, so graviti- 
nos may be much lighter than all other sparticles, with implications for 
both collider physics and cosmology. In anomaly-mediated SUSY break- 
ing (AMSB) models, 14 gravitational anomalies induce SSB terms, and the 
gravitino can be much heavier than the weak scale. There are yet other 
models 15 where SSB parameters get comparable contributions from gravity- 
mediated as well as from anomaly-mediated sources, and very recently, also 
from gauge- mediation. 16 The pattern of superpartner masses is sensitive 
to the mediation-mechanism, so that we can expect collider experiments 
to reveal which of the various mechanisms that have been proposed are 
actually realized in nature. We also mention that in both the GMSB and 
AMSB models, it is somewhat less natural (but still possible!) to obtain the 
required amount of SUSY dark matter in the Universe. Although these are 
all viable scenarios, they have not been as well scrutinized as the mSUGRA 
model. 

3. Lecture 2: Sparticle production, decay and event 
generation 

3.1. The role of the CERN Large Hadron Collider 

The CERN Large Hadron Collider (LHC) turned on briefly in 2008 before 
an electrical mis-connection caused a helium leak, and a shut-down until 
mid-2009. The LHC is a pp collider, expected to operate ultimately at 
center-of-mass energy y/s = 14 TeV. Since protons are not fundamental 
particles, we may think instead of LHC as being a quark-gluon collider, 
where the quark and gluon constituents of the proton that undergo hard 
strong or electroweak scatterings only carry typically a small fraction of the 
proton energy. 

For the hadronic reaction 

A + B^c + d + X, (10) 

where A and B are hadrons, c and d might be superpartners, and X is 
associated hadronic debris, we are really interested in the subproces reaction 

a + b^c + d, (11) 

where a is a partonic constituent of hadron A, and b is a partonic constituent 
of hadron B. The SM or MSSM or whatever effective theory we are working 
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in provides us with the Feynman rules to calculate the subprocess cross- 
section a(ab — > cd). To find the total hadron-hadron cross section, we must 
convolute the subprocess cross section with the parton distribution functions 
fa/A( x a, Q 2 ) which gives the probability to find parton a in hadron A with 
momentum fraction x a at energy-squared scale Q 2 . Thus, we have 

da(AB -» cdX) = V / dx a dx b f a/A (x a ,Q 2 ) f b/B {x b , Q 2 )da{ab -> cd). 

o,6 J ° J ° 

(12) 

where the sum extends over all initial partons a and b whose collisions 
produce the final state particles c and d. The parton-parton center-of-mass 
energy-squared s = x a x b s which enters the hard scatterng is only a small 
fraction of the proton-proton center-of-mass energy-squared s. Thus, to 
explore the TeV scale, a hadron-hadron collider of tens of TcV is indeed 
needed! 

Around the 17 mile circumference LHC ring are situated four experi- 
ments: Atlas, CMS, LHC-B and Alice. Atlas and CMS arc all-purpose de- 
tectors designed to detect almost all the energy emitted in hard scattering 
reactions. They will play a central role in the search for new physics at the 
LHC, and the search for dark matter production in the LHC environment. 

3.2. Cross section calculations for the LHC 

As noted above, if a perturbative model for new physics is available, then 
the relativistic-quantum mechanical amplitude M. for the relevant scatter- 
ing sub-process may be calculated, usually at a low order in perturbation 
theory. The probability for scattering is given by the square of the (com- 
plex) scattering amplitude M.M) , and can often be expressed as a function 
of dot products of the external leg momenta entering the diagrams, after 
summing and averaging over spins and colors. Nowadays, several computer 
programs- such as CalcHEP, CompHEP and MadGraph- are available to 
automatically do tree level calculations of various 2 — > n processes. 

To gain the parton level scattering cross section, one must multiply by 
suitable phase space factors and divide by the flux factor. The phase space 
integrals for simple sub-processes may be done analytically. More often, 
the final integrals arc done using Monte Carlo (MC) techniques, which effi- 
ciently allow integration over multi-body phase space. The MC technique, 
wherein random integration variables are generated computationally, and 
summed over, actually allows for a simulation of the scattering process. 
This can allow one to impose cuts, or simulate detectors, to gain a more 
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exact correspondence to the experimental environment. 
3.3. Sparticle production at the LHC 

Direct production of neutralino dark matter at the LHC (pp — > Z1Z1X, 
where X again stands for assorted hadronic debris) is of little interest since 
the high px final state particles all escape the detector, and there is little 
if anything to trigger an event record. Detectable events come from the 
production of the heavier supcrpartners, which in turn decay via a multi- 
step cascade which ends in the stable lightest SUSY particle (LSP). 

In many SUSY models, the strongly interacting squarks and/or gluinos 
are among the heaviest states. Unless these are extremely heavy, they will 
have large production cross sections at the LHC. Strong interaction pro- 
duction mechanisms for their production include, 1. gluino pair production 
gg, 2. squark pair production qq and 3. squark-gluino associated production 
qg. Note here that the reactions involving squarks include a huge number 
of subprocess reactions to cover the many flavors, types (left- and right- 
), and also the anti-squarks. The various possibilities each have different 
angular dependence in the production cross sections, 17 and the different 
flavors/types of squarks each have different decay modes. 18 These all have 
to be kept track of in order to obtain a reliable picture of the implications 
of SUSY in the LHC detector environment. Squarks and gluinos can also 
be produced in association with charginos and neutralinos. 19 Associated 
gluino production occurs via squark exchange in the t or u channels and is 
suppressed if squarks are very heavy 

If colored sparticles are very heavy, then electroweak production of 
charginos and neutralinos may be the dominant sparticle production mech- 
anism at the LHC. The most important processes are pair production of 
charginos, W^Wj where i,j = 1,2, and chargino-neutralino production, 
W^Zj, with « = 1,2 and j = 1 — 4. In models with unified GUT scale gaug- 
ino masses and large \u\, ZW\W\ and WZ 2 Wi couplings are large so that 
WiWi and W\Z 2 production occurs at significant rates. The latter process 
can lead to the gold-plated trilcpton signature at the LHC. 20 Neutralino 
pair production (pp — > ZiZjX where i, j = 1 — 4) is also possible. This 
reaction occurs at low rates at the LHC unless \u\ ~ M\ t 2 ( as in the case 
of mixed higgsino dark matter (MHDM)). Finally, we mention slepton pair 
production: ( + (r , v^t and vtve, which can give detectable dilepton signals 
if m- t ~ 300 GeV. 21 

In Fig. 2 we show various sparticle production cross sections at the LHC 
as a function of m§. Strong interaction production mechanisms dominate 
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at low mass, while electroweak processes dominate at high mass. The as- 
sociated production mechanisms are never dominant. The expected LHC 
integrated luminosity in the first year of running is expected to be around 
0.1 fb _1 , while several tens of fb _1 of data is expected to be recorded in the 
first several years of operation. The ultimate goal is to accumulate around 
500-1000 fb-\ correponding to 10 5 - 10 6 SUSY events for m g ~ 1 TeV. 
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Fig. 2. Cross sections for production of various sparticles at the LHC. Gaugino 
mass unification is assumed. 



3.4. Sparticle cascade decays 

In imparity conserving SUSY models, sparticles decay to lighter sparticles 
until the decay terminates in the LSP. 18 Frequently, the direct decay to 
the LSP is either forbidden or occurs with only a small branching fraction. 
Since gravitational interactions are negligible, gluinos can only decay via 
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g — > qq, where the q and q can be of any flavor or type. If two body decay 
modes are closed, the squark will be virtual, and the gluino will decay via 
three body modes g — > qqZ i7 qq'Wj. If squarks are degenerate, and Yukawa 
coupling effects negligible, three-body decays to the wino-like chargino and 
neutralino usually have larger branching fractions on account of the larger 
gauge coupling. If < M 2 , gluinos and squarks may thus decay most 
of the time to the heavier charginos and neutralinos, resulting in lengthy 
cascade decay chains at the LHC. 

Squarks decay always to two-body modes: q — > qg if it is kincmatically 
allowed, or q~L — > q'W, qZj, while q~R — > qZj only, since right-squarks do 
not couple to charginos. Sleptons do not have strong interactions so cannot 
decay to gluinos. Their electroweak decays are similar to corresponding 
decays of squarks II — * i'Wi, IZj while Ir — ► £Zj only. 

Charginos may decay via two-body modes: Wi — > WZ,-, ZWj 
or even to </>Wj or H~Zj, where <j> — h,H,A. If two-body modes are in- 
accessible, then three-body decays dominate: Wi — > Zjff', where / and 
/' are SM fermions which couple to the W. Frequently, the decay ampli- 
tude is dominated by the virtual W so that the three-body decays of W\ 
have the same branching fractions as those of the W. Neutralinos decay via 
Zi — > WWj, H + Wj, ZZj, <f>Zj or //. If two body neutralino decays are 
closed, then Zi — > Zjff, where / are the SM fermions. In some models, 
the branching fraction for radiative decays Zi — ► Zjj (that only occurs 
at the one-loop level) may be significant. 22 The cascade decay modes of 
neutralinos depend sensitively on model parameters. 23 

If tan (3 is large, then b and r Yukawa coupling effects become important, 
enhancing three body decays of g, Wi and Zj to third generation fermions. 24 
For very large values of tan j3 these decays can even dominate, resulting in 
large rates for 6-jet and r-jet production in SUSY events. 25 

Finally, the various Higgs bosons can be produced both directly and via 
sparticle cascades at the LHC. 26 Indeed, it may be possible that h is first 
discovered in SUSY events because in a sample of events enriched for SUSY, 
it is possible to identify h via its dominant h — > bb decays rather than via 
its sub-dominant decay modes, as required for conventional searches. 26 The 
heavier Higgs bosons decay to a variety of SM modes, but also to SUSY 
particles if these latter decays are kincmatically allowed, leading to novel 
signatures such a.s H, A ^> Z 2 Z 2 — ► 4^+ $t- 27 

The cascade decays terminate in the LSP. In the case of a Z\ LSP, the 
Z\ is a DM candidate, and leaves its imprint via $t- In the case of a weak 
scale G or a LSP, then Z\ will decay as Z\ — ► 7G or 70. In these cases, 
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the Z\ lifetime is long enough that it decays outside the detector, so one 
still expects large fix in the collider events. An exception arises for the 
case of super-light gravitinos (with masses in the eV to keV range) that are 
possible in GMSB models: see (15). Then, the decay may take place inside 
the detector, possibly with a large vertex separation. It is also possible that 
the NLSP is charged and quasi-stable, in which case collider events may 
include highly ionizing tracks instead of, or in addition to, Et- 

The decay branching fractions depend on the entire spectrum of SUSY 
particle masses and and their mixings. They are pre-programmed in several 
codes: Isajet, 28 SDECAY 29 and Spheno. 30 

3.5. Event generation for LHC 

Once sparticle production cross sections and decay branching fractions have 
been computed, it is useful to embed these into event generator programs 
to simulate what SUSY collider events will look like at LHC. There are 
several steps involved: 

• Calculate all sparticle pair production cross sections. Once all ini- 
tial and final states are accounted for, this involves over a thousand 
individual subprocess reactions. In event generation, a particular 
reaction is selected on a probabilistic basis, with a weight propor- 
tional to its differential cross-section. 

• Sparticle decays are selected probabilistically into all the allowed 
modes in proportion to the corresponding branching fractions. 

• Initial and final state quark and gluon radiation are usually dealt 
with using the parton shower (PS) algorithm, which allows for prob- 
abilistic parton emission based on approximate collincar QCD emis- 
sion matrix elements, but exact kinematics. The PS is also applied 
at each step of the cascade decays, which may lead to additional 
jet production in SUSY collider events. 

• A hadronization algorithm provides a model for turning various 
quarks and gluons into mesons and baryons. Unstable hadrons must 
be further decayed. 

• The beam remnants - proton constituents not taking part in the 
hard scattering - must be showered and hadronized, usually with 
an independent algorithm, so that energy deposition in the forward 
detector region may be reliably calculated. 

At this stage, the output of an event generator program is a listing of 
particle types and their associated four-vectors. The resulting event can 
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then be interfaced with detector simulation programs to model what the 
actual events containing DM will look like in the environment of a collider 
detector. 

Several programs are available, including Isajet, 28 Pythia, 31 Herwig 32 
and Sherpa. 33 Other programs such as Madevent, 34 CompHEP/CalcHEP 35 
and Whizard 36 can generate various 2 — > n processes including SUSY par- 
ticles. The output of these programs may then be used as input to Pythia 
or Herwig for showering and hadronization. Likewise, parton level Isajet 
SUSY production followed by cascade decays can be input to Pythia and 
Herwig via the Les Houches Event format. 37 

4. Lecture 3: SUSY, LHT and UED signatures at the LHC 
4.1. Signatures for SUSY particle production 

Unless colored sparticles are very heavy, the SUSY events at the LHC 
mainly result in gluino and squark production, followed by their possi- 
bly lengthy cascade decays. These events, therefore, typically contain very 
hard jets (from the primary decay of the squark and/or gluino) together 
with other jets and isolated electrons, muons and taus (identified as narrow 
one- and three-prong jets), and sometimes also photons, from the decays 
of secondary charginos and neutralinos, along with fix that arises from 
the escaping dark matter particles (as well as from neutrinos). In models 
with a superlight gravitino, there may also be additional isolated photons, 
leptons or jets from the decay of the NLSP. The relative rates for vari- 
ous n-jet + m-lcpton + fc-photon +$t event topologies is sensitive to the 
model as well as to the parameter values, and so provide a useful handle 
for phenomenological analyses. 

Within the SM, the physics background to the classic jets+ fix sig- 
nal comes from neutrinos escaping the detector. Thus, the dominant SM 
backgrounds come from W + jets and Z + jets production, ti production, 
QCD multijet production (including bb and cc production), WW, WZ, ZZ 
production plus a variety of 2 — > n processes which are not usually included 
in event generators. These latter would include processes such as titi, ttbb, 
tiW, WWW, WWZ production, etc. Decays of electroweak gauge bosons 
and the i-quark are the main source of isolated leptons in the SM. Var- 
ious additional effects- uninstrumented regions, energy mis- measurement, 
cosmic rays, beam-gas events- can also lead to $t events. 

In contrast to the SM, SUSY events naturally tend to have large jet 
multiplicities and frequently an observable rate for high multiplicity lepton 
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events with large J/?t- Thus, if one plots signal and background versus mul- 
tiplicity of any of these quantities, as one steps out to large multiplicity, 
the expected SUSY events should increase in importance, and even dom- 
inate the high multiplicity channels in some cases. This is especially true 
of isolated multi-lepton signatures, and in fact it is convenient to classify 
SUSY signal according to lepton multiplicity: 38 

• zero lepton +jets+ $t events, 

• one lepton +jets+ $t events, 

• two opposite sign leptons +jets+$r events (OS), 

— same- flavor (OSSF), 

- different flavor (OSDF), 

• two same sign leptons +jets+ $t events (SS), 

• three leptons +jets+ $t events (31), 

• four (or more) leptons +jets+ $t events (4£). 

4.2. LHC reach for SUSY 

Event generators, together with detector simulation programs can be used 
to project the SUSY discovery reach of the LHC. Given a specific model, one 
may first generate a grid of points that samples the parameter (sub)space 
where signal rates are expected to vary significantly. A large number of 
SUSY collider events can then be generated at every point on the grid along 
with the various SM backgrounds to the SUSY signal mentioned above. 
Next, these signal and background events are passed through a detector 
simulation program and a jet- finding algorithm is implemented to determine 
the number of jets per event above some Er(jei) threshold (usually taken 
to be Er(jei) > 50 — I 00 GeV for LHC). Finally, analysis cuts are imposed 
which are designed to reject mainly SM BG while retaining the signal. 
These cuts may include both topological and kinematic selection criteria. 
For observability with an assumed integrated luminosity, we require that 
the signal exceed the chance 5 standard deviation upward fluctuation of 
the background, together with a minimum value of (~ 25%) the signal to 
background ratio, to allow for the fact that the background is not perfectly 
known. For lower sparticle masses, softer kinematic cuts are used, but for 
high sparticle masses, the lower cross sections but higher energy release 
demand hard cuts to optimize signal over background. 

In Fig. 3, we illustrate the SUSY reach of the LHC within the mSUGRA 
model assuming an integrated luminosity of 100 fb _1 . We show the result 
in the toq — plane, taking A = 0, tan/3 = 10 and \i > 0. The signal 
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is observable over background in the corresponding topology below the 
corresponding curve. We note the following. 

(1) Unless sparticles are very heavy, there is an observable signal in several 
different event topologies. This will help add confidence that one is 
actually seeing new physics, and may help to sort out the production 
and decay mechanisms. 

(2) The reach at low mo extends to m 1 / 2 ~ 1400 GeV. This corresponds 
to a reach for mq ~ rrig ^3.1 TeV. 

(3) At large mo, squarks and sleptons are in the 4—5 TeV range, and are 
too heavy to be produced at significant rates at LHC. Here, the reach 
comes mainly from just gluino pair production. In this range, the LHC 
reach is up to m 1( / 2 ~ 700 GeV, corresponding to a reach in m§ of 
about 1.8 TeV, and may be extended by ~ 15-20% by 6-jet tagging. 39 

mSuera with tanPt = 30, A~ = 0, ll > 
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Fig. 3. The 100 ftT 1 fb reach of LHC for SUSY in the mSUGRA model. For 
each event topology, the signal is observable below the corresponding contour. 

In Fig. 4 we can see a comparison of the LHC reach (notice that it is 
insensitive to tan/3 and sign(/i)) with that of the Tevatron (for clean 2,1 
events with 10 fb -1 ), and the proposed e + e~ International Linear Collider 
(ILC), with ,/s = 0.5 or 1 TeV along with various dark matter direct 
detection (DD) and indirect detection (ID) search experiments. We remark 
that: 
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• While LHC can cover most of the relic density allowed region, the 
HB/FP region emerges far beyond the LHC reach. 

• As will be discussed, the DD and ID experiments have the greatest 
sensitivity in the HB/FP region where the neutralino is MHDM. In this 
sense, DD and ID experiments complement LHC searches for SUSY. 

• The ILC reach is everywhere lower than LHC, except in the HB/FP re- 
gion. In this region, while gluinos and squarks can be extremely heavy, 
the a parameter is small, leading to a relatively light spectrum of 
charginos and neutralinos. These are not detectable at the LHC because 
the visible decay products are too soft. However, since chargino pair 
production is detectable at ILC even if the energy release in chargino 
decays is small, the ILC reach extends beyond LHC in this region. 40 

Finally, we note here that while the results presented above are for the 
LHC reach in the mSUGRA model, the LHC reach (measured in terms of 
wig and niq) tends to be relatively insensitive to the details of the model 
chosen, as long as gluino and squark production followed by cascade decays 
to the DM particle occur. 

4.3. Early discovery of SUSY at LHC without Et 

The classic signature for SUSY at hadron colliders is the presence of jets, 
isolated leptons and especially large Et- However, in the early period of 
LHC running, fix may be very difficult to measure reliably. In a real de- 
tector, there will be calorimetry calibration issues, dead or uninstrumented 
regions, "hot" cells, cosmic ray events, beam-gas collisions: all of these con- 
tribute to Et events, in additon to SM backgrounds and possibly new 
physics. As learned from Tevatron experiments, a reliable measurement of 
Et requires knowledge of the complete detector, and that may well take 
some time- possibly over a year- after start-up. 

Can LHC search for SUSY even if Et searches are not viable? The an- 
swer seems to be yes, as long as lepton isolation cuts are possible. In Rcf., 41 
it was shown that requiring events with > 4 jets plus a high multiplicity 
of isolated leptons: SS, OSSF, OSDF, 3£, ■ ■ ■ would allow one to severely 
reduce SM background while maintaining large enough signal rates. It is 
claimed that an LHC reach of m g ~ 750 GeV is possible with just 0.1 fb _1 
of integrated luminosity, without using an $t cut, by requiring events with 
> 3 isolated leptons (es or as). 

In addition, reliable electron identification may also be an issue in early 
LHC running. The problem here is differentiating the electron's EM shower 



January 29, 2009 20:12 



WSPC - Proceedings Trim Size: 9in x 6in tasi 



21 

mSUGRA, A =0, tan|3=45, n<0 

1600 
1400 
1200 
§1000 
E 800 
600 
400 
200 

1000 2000 ,„ ,,,3000 4000 5000 

rm (GeV) 

4>(p")=3e" 7 GeV 1 cm" 2 s" 1 sr"' _ (S/B) et =0.01 
_4(Y)=10" 1 °cm" 2 s" 1 _ 4> Bun (n)=40 km" 2 yr"' m h =1 14.4 GeV 
_ <5""'"'((i)=40 km" 2 yr"' _ a{Zfi)=10 ! ' pb 

0< !lh 2 < 0.129 

Fig. 4. The projected reach of various colliders, direct and indirect dark matter 
search experiments in the mSUGRA model. For the indirect search results we have 
adopted the conservative default DarkSUSY isotropic DM halo density distribu- 
tion. Plot is from Ref. 87 



from other energy depositions such as charged pions, or hard photons where 
in additon a charged track is nearby. In this case, a SUSY search may 
still be made by looking for multi-jet plus multi-muon production, again 
without any requirement on missing Et- A2 Atlas and CMS are already 
reliably detecting cosmic muons. Also, muons enjoy an addvantage over 
electrons in that they can be detected reliably to pr values as low as 5 
GeV. 

The multiplicity of muons in > 4 jet events is shown in Fig. 5, assuming 
Pt(m) > 10 GeV to remove BG from heavy flavor decays in SM processes. 
Here, we see that as one moves to high muon multiplicity, the SM back- 
ground drops off sharply, and by > 3, signal already far exceeds BG, 
at least for the SPSla' signal point shown in the figure. In fact, by requir- 
ing same-sign dimuon plus multi-jet events, already signal in many cases 
will exceed background. The LHC should be able to discover gluinos with 
rrig ~ 450 (550) GeV in the SS dimuon plus > 4 jets state with just 0.04 
(0.1) fb _1 of integrated luminosity. 
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Fig. 5. Multiplicity of isolated muons in> 4 jets events, from signal point SPSla' 
and various SM background processes. 



4.4. Determination of sparticle properties 

Once a putative signal for new physics emerges at LHC, the next step is to 
establish its origin. This will entail detailed measurements of cross sections 
and distributions in various event topologies to gain insight into the identity 
of the new particles being produced, their masses, decay patterns, spins, 
couplings (gauge quantum numbers) and ultimately mixing angles. These 
measurements are not straightforward in the LHC environment because 
of numerous possible SUSY production reactions occurring simultaneously, 
a plethora of sparticle cascade decay possibilities, hadronic debris from 
initial state radiation and lack of invariant mass reconstruction due to the 
presence of fir- All these lead to ambiguities and combinatoric problems 
in reconstructing exactly what sort of signal reactions are taking place. In 
contrast, at the ILC, the initial state is simple, the beam energy is tunable 
and beam polarization can be used to select out specific processes. 

While it seems clear that the ILC is better suited for a systematic pro- 
gram of precision sparticle measurements, studies have shown (albeit in 
special cases) that interesting measurements are also possible at the LHC. 
We go into just a subset of all details here in order to give the reader an 
idea of some of the possibilities suggested in the literature. 

One suggested starting point is the distribution of effective mass M c g = / 
Erijets) in the inclusive SUSY sample, which sets the approximate 
mass scale Msusy = min(mg, rriq) for the strongly interacting sparticles are 
being produced, 43 and provides a measure of Msusy to 10-15%. 
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More detailed information on sparticle masses may be accessed by 
studying specific event topologies. For instance, the mass of dileptons from 
Zi — > £ + £r Z\ decays is bounded by — (this bound is even more 

restrictive if Zi decays via an on-shell slepton). 44 We therefore expect an 
OSSF invariant mass distribution to exhibit an edge at m^ 2 — m z 1 ( or ^° c ~ 
low) in any sample of SUSY events so long as the "spoiler" decay modes 
Z2 — * Z\Z or Z\h are closed. Contamination from chargino production can 
be statistically removed by subtracting out the distribution of OSDF dilep- 
tons. In MHDM models, there may be more than one visible mass edge 
because the Z% may also be accessible in cascade decays. 

In the happy circumstance where production of gluinos or a single type 
of squark is dominant, followed by a string of two-body decays, then further 
invariant mass edges are possible. One well-studied example comes from 
g — > bbi — > bbZ2 — >• bb££Z\, then one can try to combine a 6-jet with 
the dilepton pair to reconstruct the squark-neutralino mass edge: m(b££) < 
m b 1 ~ m z L - Next, combining with another 6-jet can yield a gluino-neutralino 
edge: m(bb££) < trig — . The reconstruction of such a decay chain may be 
possible as shown in Ref., 43 where other sequences of two-body decays are 
also examined. In practice, such fortuitous circumstances may not exist, and 
there are many combinatoric issues to overcome as well. A different study 45 
shows that end-point measurements at the LHC will make it possible to 
access the mass difference between the LSP and the stau in a mSUGRA 
scenario where the stau co-annihilation mechanism is operative. 

These end-point measurements generally give mass differences, not 
masses. However, by an analysis of the decay chain q~L — > <Z-^2 — » qi^t^ — > 
qi^t^Zi, it has been argued 46 that reconstruction of masses may be possi- 
ble under fortuituous circumstances. More recently, it has been suggested 
that it may be possible to directly access the gluino and/or squark masses 
(not mass differences) via the introduction of the so-called niT2 variable. 
We will refer the reader to the literature for details. 47 It also may be pos- 
sible to make a measurement of m g basd on total production cross section 
in cases where pure gluino pair production is dominant. 48 

Mass measurements allow us to check consistency of specific SUSY mod- 
els with a handful of parameters, and together with other measurements 
can readily exclude such models. But these are not the only interesting 
measurements at the LHC. It has been shown that if the NLSP of GMSB 
models decays into a superlight gravitino, it may be possible to determine 
its lifetime, and hence the gravitino mass at the LHC. 49 This will then 
allow one to infer the underlying SUSY breaking scale, a scale at least as 
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important as the weak scale! A recent study 50 suggests that this is possible 
even when the the decay length of the NLSP is too short to be measured. 
While linear collider experiments will ultimately allow the precision mea- 
surements that will directly determine the new physics to be softly broken 
supersymmctry, 51 it will be exciting to analyze real LHC data that will 
soon be available to unravel many of the specific details about how (or if) 
SUSY is actually implemented in nature. 

4.5. Measuring DM properties at LHC and ILC 

SUSY discovery will undoubtedly be followed by a program (as outlined 
in Sec. 4.4) to reconstruct sparticle properties. What will we be able to 
say about dark matter in light of these measurements? Such a study was 
made by Baltz et al. 52 where four mSUGRA case study points (one each in 
the bulk region, the HB/FP region, the stau coanihilation region and the 
^4-funnel region) were examined for the precision with which measurements 
of sparticle properties that could be made at LHC, and also at a ^/s = 0.5 
and 1 TeV e + e~ collider. They then adopted a 24-parameter version of 
the MSSM and fit its parameters to these projected measurements. The 
model was then used to predict several quantities relevant to astrophysics 
and cosmology: the dark matter relic density fir, the spin- independent 
neutralino-nucleon scattering cross section asi{Z\p), and the neutralino 
annihilation cross section times relative velocity, in the limit that v — > 0: 
(av)\ v —>o- The last quantity is the crucial particle physics input for estimat- 
ing signal strength from neutralino annihilation to anti-matter or gammas 
in the galactic halo. What this yields then is a collider measurement of 
these key dark matter quantities. 

As an illustration, we show in Fig. 6 (taken from Ref. 52 ) the precision 
with which the neutralino relic density is constrained by collider measure- 
ments for the LCC2 point which is in the HB/FP region of the mSUGRA 
model. Measurements at the LHC cannot fix the LSP composition, and so 
are unable to resolve the degeneracy between a wino-LSP solution (which 
gives a tiny relic density) and the true solution with MHDM. Determina- 
tions of chargino production cross sections at the ILC can easily resolve the 
difference. It is nonetheless striking that up to this degeneracy ambiguity, 
experiments at the LHC can pin down the relic density to within ~ 50% (a 
remarkable result, given that there are sensible models where the predicted 
relic density may differ by orders of magnitude!). This improves to 10-20% 
if we can combine LHC and ILC measurements. 
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This collider determination of the relic density is very important. If it 
agrees with the cosmological measurement it would establish that the DM 
is dominantly thermal ncutralinos from the Big Bang. If the neutralino 
relic density from colliders falls significantly below (1), it would provide 
direct evidence for multi-component DM- perhaps neutralinos plus axions 
or other exotica. Alternatively, if the collider determination gives a much 
larger value of f2 ^ h 2 , it could point to a long-lived but unstable neutralino 
and/or non-thermal DM. 

The collider determination of model parameters would also pin down 
the neutralino-nucleon scattering cross section. Then if a WIMP signal is 
actually observed in DD experiments, one might be able to determine the 
local DM density of neutralinos and aspects of their velocity distribution 
based on the DD signal rate. This density should agree with that obtained 
from astrophysics if the DM in our Galaxy is comprised only of neutralinos. 

Finally, a collider determination of (av)\ v —to would eliminate uncer- 
tainty on the particle physics side of projections for any ID signal from 
annihilation of neutralinos in the galactic halo. Thus, the observation of 
a gamma ray and/or anti-matter signal from neutralino halo annihilations 
would facilitate the determination of the galactic halo dark matter density 
distribution. 
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4.6. Some non-SUSY WIMPs at the LHC 

4.6.1. B^ state from universal extra dimensions 

Models with Universal Extra Dimensions, or UED, are interesting con- 
structs which provide a foil for SUSY search analyses. 53 In the 5-D UED 
theory, one posits that the fields of the SM actually live in a 5-D brane 
world. The extra dimension is "universal" since all the SM particles prop- 
agate in the 5-D bulk. The single extra dimension is assumed to be com- 
pactified on a S1/Z2 orbifold (line segment). After compactification, the 
4-D effective theory includes the usual SM particles, together with an infi- 
nite tower of Kaluza-Klein (KK) excitations. The masses of the excitations 
depend on the radius of the compactified dimension, and the first (n = 1) 
KK excitations can be taken to be of order the weak scale. In these theories, 
KK-parity (—1)™ can be a conserved quantum number. If this so-called KK- 
parity is exact, then the lightest odd KK parity state will be stable and can 
be a DM candidate. At tree-level, all the KK excitations in a given level are 
essentially degenerate. Radiative corrections break the degeneracy, leaving 
colored excitations as the heaviest excited states and the n = 1 KK excita- 
tion of the SM U(l)y gauge boson as the lightest 54 KK odd state: in the 
UED case, therefore, the DM particle has spin-1. The splitting caused by 
the radiative corrections is also essential to assess how the KK excitations 
decay, and hence are crucial for collider phenomenology. 55 

The relic density of B* particles has been computed, and found to be 
compatible with observation for certain mass ranges of B^. 56 Also, in UED, 
the colored excitations can be produced with large cross sections at the 
LHC, and decay via a cascade to the final state. Thus, the collider 
signatures are somewhat reminiscent of SUSY, and it is interesting to ask 
whether it is possible to distinguish a jets+leptons+ $t signal in UED from 
that in SUSY. Several studies 57 answer affirmatively, and in fact provide 
strong motivation for the measurement of the spins of the produced new 
particles. 58 UED DM generally leads to a large rate in IceCube, and may 
also give an observable signal in anti-protons and possibly also in photons 
and positrons. 53,59 DD is also possible but the SI cross section is typically 
smaller than 1CU 9 pb. 

4.6.2. Little Higgs models 

Little Higgs models 60,61 provide an alternative method compared to SUSY 
to evade the quadratic sensitivity of the scalar Higgs sector to ultra-violet 
(UV) physics. In this framework, the Higgs boson is a pseudo-Goldstonc 
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boson of a spontaneously broken global symmetry that is not completely 
broken by any one coupling, but is broken when all couplings are included. 
This then implies that there is quadratic sensitivity to UV physics, but 
only at the multi-loop level. Specific models where the quadratic sensitivity 
enters at the two-loop level should, therefore, be regarded as low energy 
effective theories valid up to a scale A ~ 10 TeV, at which a currently 
unknown, and perhaps strongly-coupled UV completion of the theory is as- 
sumed to exist. Models that realize this idea require new TeV-scale degrees 
of freedom that can be searched for at the LHC: new gauge bosons, a heavy 
top-like quark, and new spin-zero particles, all with couplings to the SM. 
These models, however, run into phenomenological difficulties with preci- 
sion EW constraints, unless a discrete symmetry- dubbed T-parity 62 - is 
included. SM particles are then T-even, while the new particles are T-odd. 

We will set aside the issue of whether T-parity conservation is violated 
by anomalies, 64 and assume that a conserved T-parity can be introduced. 65 
In this case, the lightest T-odd particle Ah - the Little Higgs partner of 
the hypercharge gauge boson with a small admixture of the neutral W^h 
boson - is stable and yields the observed amount of DM for a reasonable 
range of model parameters. 59 In this case, the DM particle has spin-1, 
though other cases with either a spin-i or spin-0 heavy particle may also 
be possible. Ah can either annihilate with itself into vector boson pairs or 
tt pairs via s-channel Higgs exchange, or into top pairs via exchange of the 
heavy T-odd quark in the t-channel. Co-annihilation may also be possible 
if the heavy quark and Ah are sufficiently close in mass. Signals at the 
LHC 66 ' 67 mainly come from pair production of heavy quarks pp — > TTX 
followed by T — > tAn decay, and from single production of the heavy quark 
in association with Ah- These lead to tt+ fix events at LHC. 68 The fix 
comes from the escaping Ah particle, which must be the endpoint of all 
T-odd particle decays. a If Ah is the dominant component of galactic DM, 
we will generally expect small DD and ID rates for much the same reasons 
that the signals from the bino LSP tend to be small: 59 see, however, Ref. 70 
for a different model with large direct detection rate. 



a We note here that it is also possible to construct so-called twin-Higgs models 69 where 
the Higgs sector is stabilized via new particles that couple to the SM Higgs doublet, 
but are singlets under the SM gauge group. In this case, there would be no obvious new 
physics signals at the LHC. 
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5. Lecture 4: $t and the dark matter connection 
5.1. Neutralino relic density 

Once a SUSY model is specified, then given a set of input parameters, it 
is possible to compute all superpartner masses and couplings necessary for 
phenomenology. We can then use these to calculate scattering cross sections 
and sparticle decay patterns to evaluate SUSY signals (and corresponding 
SM backgrounds) in collider experiments. We can also check whether the 
model is allowed or excluded by experimental constraints, either from direct 
SUSY searches, e.g. at LEP2 which requires that > 103.5 GeV, mg ~ 
100 GeV, and m h > 114.4 GeV (for a SM-like light SUSY Higgs boson h), 
or from indirect searches through loop effects from SUSY particles in low 
energy measurements such as B(b — + sj) or (g — 2) M . We can also calculate 
the expected thermal LSP relic density. To begin our discussion, we will 
first assume that the lightest neutralino Z\ is the candidate DM particle. 

As mentioned earlier, the relic density calculation involves solving the 
Boltzmann equation, where the neutralino density changes due to both the 
expansion of the Universe and because of neutralino annihilation into SM 
particles, determined by the thermally averaged Z\Z\ annihilation cross 
section. An added complication occurs if neutralino co- annihilation is pos- 
sible. Co-annihilation occurs if there is another SUSY particle close in mass 
to the Z\ , whose thermal relic density (usually suppressed by the Boltzmann 
factor exp ~^ M ) is also significant. In the mSUGRA model, co- annihilation 
may occur from a stau, fi , a stop t\ or the lighter chargino W\ . For instance, 
in some mSUGRA parameter-space regions the f\ and Z\ are almost degen- 
erate, so that they both have a significant density in the early universe, and 
reactions such as Z\f\ — ► occur. Since the electrically charged fi can 
also annihilate efficiently via electromagnetic interactions, this process also 
alters the equilibrium density of neutralinos. All in all, there are well over 
a thousand neutralino annihilation and co-annihilation reactions that need 
to be computed, involving of order 7000 Feynman diagrams. There exist 
several publicly available computer codes that compute the neutralino relic 
density: these include DarkSUSY, 71 MicroMegas 72 and IsaRcD 73 (a part of 
the Isatools package of Isajet 28 ). 

As an example, we show in Fig. 7 the m n vs. mi/2 plane from the 
mSUGRA model, where we take A a = 0, fi > 0, m t — 171.4 GeV and 
tan/3 = 10. The red-shaded regions arc not allowed because either the fi 
becomes the lightest SUSY particle, in contradiction to negative searches 
for long lived, charged relics (left edge), or EWSB is not correctly obtained 
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(lower-right region). The blue-shaded region is excluded by LEP2 searches 
for chargino pair production [rn^ < 103.5 GeV). We show contours of 
squark (solid) and gluino (dashed) mass (which are nearly invariant under 
change of A$ and tan 0) . Below the magenta contour near miw ~ 200 GeV, 
rrih < 110 GeV, which is roughly the LEP2 lower limit on mu in the model. 
The thin green regions at the edge of the unshaded white region has 0.094 < 
fijj h 2 < 0.129 where the neutralino saturates the observed relic density. 
In the adjoining yellow regions, 0.% h 2 < 0.094, so these regions require 
multiple DM components. The white regions all have il^ h 2 > 0.129 and 
so give too much thermal DM: they are excluded in the standard Big Bang 
cosmology. 



mSUGRA : tan|3=10, A =0, n>0, m,=171.4 GeV 




1 2 3 4 5 



m (TeV) 

0<ilh 2 < 0.094 • Excluded 

• 0.094 <£2h 2 < 0.129 • LEP2 

Fig. 7. DM-allowed regions in the tuq — m 1 y 2 plane of the mSUGRA model for 
tan p — 10 with Aq = and yi > 0. 

The DM-allowed regions are classified as follows: 

• At very low mo and low m-yn values is the so-called bulk annihilation 
region. 74 Here, sleptons are quite light, so Z\Zi ~ ■> it via t-channcl 
slepton exchange. In years past (when ^cdm^ 2 ~ 0.3 was quite con- 
sistent with data), this was regarded as the favored region. But today 
LEP2 sparticle search limits have increased the LEP2-forbidden region 
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from below, while the stringent bound flcDMh 2 < 0.13 has pushed the 
DM-allowed region down. Now hardly any bulk region survives in the 
mSUGRA model. 

• At low m and moderate r/iiw, there is a thin strip of (barely dis- 
cernable) allowed region adjacent to the stau-LSP region where the 
neutralino and the lighter stau were in thermal equilibrium in the early 
universe. Here co-annihilation with the light stau serves to bring the 
neutralino relic density down to its observed value. 75 

• At large m , adjacent to the EWSB excluded region on the right, is 
the hyperbolic branch/focus point (HB/FP) region, where the super- 
potential fj, parameter becomes small and the higgsino-content of Z\ 
increases significantly. Then Z\ can annihilate efficiently via gauge cou- 
pling to its higgsino component and becomes mixed higgsino-bino DM. 
If > Mw, Mz, then ZiZi — ► WW, ZZ, Zh is enhanced, and one 
finds the correct measured relic density. 76 

We show the corresponding situation for tan/3 = 52 in Fig. 8. While the 
stau co-annihilation and the HB/FP regions are clearly visible, we see that 
a large DM consistent region now appears. 

• In this region, the value of tua is small enough so that Z\Zi can anni- 
hilate into bb pairs through s-channcl A (and also H) resonance. This 
region has been dubbed the yl-funnel. 77 It can be very broad at large 
tan (3 because the width r a can be quite wide due to the very large 
b- and r- Yukawa couplings. If tan/3 is increased further, then Z\Z\ 
annihilation through the (virtual) A* is large all over parameter space, 
and most of the theoretically-allowed parameter space becomes DM- 
consistent. For even higher tan/3 values, the parameter space collapses 
due to a lack of appropriate EWSB. 

It is also possible at low rai/2 values that a light Higgs h resonance an- 
nihilation region can occur just above the LEP2 excluded region. 78 Finally, 
if A$ is large and negative, then the t\ can become light, and <~ m^, 
so that stop-neutralino co-annihilation 79 can occur. 

Up to now, we have confined our discussion to the mSUGRA framework 
in which compatibility with (1) is obtained only over selected portions of the 
mo — mi/2 plane. The reader may well wonder what happens if we relax 
the untested universality assumptions that underlie mSUGRA. Without 
going into details, we only mention here that in many simple one-parameter 
extensions of mSUGRA where the universality of mass parameters is relaxed 
in any one of the matter scalar, the Higgs scalar, or the gaugino sectors, 
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mSUGRA : A„ = 0, u > 0, tan|3 = 52, m, = 171.4 GeV 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



m„ (TeV) 

< £2h 2 < 0.094 • LEP2 Excl. : m w , < 103.5 GeV 

• 0.094 < ffli 2 < 0.129 — m h =110GeV 

Fig. 8. DM-allowed regions in the mo — fni/2 plane of the mSUGRA model for 
tan j3 — 52 with Aq — and u > 0. The various colors of shading is as in Fig. 7. 

all points in the mo — mj/j plane become compatible with the relic density 
constraint due to a variety of mechanisms: these are catalogued in Ref. 80 
Implications of the relic density measurement for collider searches must 
thus be drawn with care. 

5.2. Neutralino direct detection 

Fits to galactic rotation curves imply a local relic density of pcdm ~ 0.3 
GeV/cm 3 . For a 100 GeV WIMP, this translates to about one WIMP per 
coffee mug volume at our location in the galaxy. The goal of DD experiments 
is to detect the very rare WIMP-nucleus collisions that should be occuring 
as the earth, together with the WIMP detector, moves through the DM 
halo. 

DD experiments are usually located deep underground to shield the 
experimental apparatus from background due to cosmic rays and ambient 
radiation from the environment or from radioactivity induced by cosmic ray 
exposure. One technique is to use cryogenic crystals cooled to near absolute 
zero, and look for phonon and ionization signals from nuclei recoiling from a 
WIMP collision. In the case of the CDMS experiment 81 at the Soudan iron 
mine, target materials include germanium and silicon. Another technique 
uses noble gases cooled to a liquid state as the target. Here, the signal is 
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scintillation light picked up by photomultiplier tubes and ionization. Tar- 
get materials include xenon, 82 argon and perhaps neon. These noble liquid 
detectors can be scaled up to large volumes at relatively low cost. They 
have the advantage of fiducialization, wherein the outer layers of the de- 
tector act as an active veto against cosmic rays or neutrons coming from 
phototubes or detector walls: only single scatters from the inner fiducial 
volume qualify as signal events. A third technique, typified by the COUPP 
experiment, 83 involves use of superheated liquids such as CF 3 I located in a 
transparent vessel. The nuclear recoil from a WIMP-nucleon collision then 
serves as a nucleation site, so that a bubble forms. The vessel is monitored 
visually by cameras. Background events are typically located close to the 
vessel wall, while neutron interactions are likely to cause several bubbles to 
form, instead of just one, as in a WIMP collision. This technique allows for 
the use of various target liquids, including those containing elements such 
as fluorine, which is sensitive to spin- dependent interactions. 

The cross section for WIMP-nucleon collisions can be calculated, and 
in the low velocity limit separates into a coherent spin-independent compo- 
nent (from scattering mediated by scalar quarks and scalar Higgs bosons) 
which scales as nuclear mass squared, and a spin-dependent component 
from scattering mediated by the Z boson or by squarks, which depends 
on the WIMP and nuclear spins. 9 The scattering cross section per nucleon 
versus mwiM p serves as a figure of merit and facilitates the comparison of 
the sensitivity of various experiments using different target materials. 

In Fig. 9, we show the spin-independent Z\p cross section versus 
for a large number of SUSY models (including mSUGRA). Every color 
represents a different model. For each model, parameters are chosen so 
that current collider constraints on sparticle masses are satisfied, and fur- 
ther, that the lightest neutralino (assumed to be the LSP) saturates the 
observed relic abundance of CDM. Also shown is the sensitivity of cur- 
rent experiments together with projected sensitivity of proposed searches 
at superCDMS, Xenon-100, LUX, WARP and at a ton-sized noble liquid 
detector. The details of the various models are unimportant for our present 
purpose. The key thing to note is that while the various models have a 
branch where a^iipZi) falls off with m^, there is another branch where 
this cross-section asymptotes to just under 10~ 8 pb. 80 > 84 > 85 Points in this- 
branch (which includes the HB/FP region of mSUGRA), are consistent with 
(1) because Z\ has a significant higgsino component. Ncutralinos with an 
enhanced higgsino content can annihilate efficiently in the early universe 
via gauge interactions. Moreover, since the spin-independent DD ampli- 
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Fig. 9. The spin-independent neutralino-proton scattering cross-section vs 

in a variety of SUSY models, compatible with collider constraints where thermally 

produced Big Bang neutralinos saturate the observed dark matter density. 



tude is mostly determined by the Higgs boson-higgsino-gaugino coupling, 
it is large in models with MHDM which has both gaugino and higgsino 
components. Thus the enhanced higgsino component of MHDM increases 
both the neutralino annihilation in the early universe as well as the spin- 
independent DD rate. The exciting thing is that the experiments currently 
being deployed- such as Xenon-100, LUX and WARP- will have the sen- 
sitivity to probe this class of models. To go further will require ton-size or 
greater target material. 

We note here that if mwiMP ~ 150 GeV, then it may be possible to 
extract the WIMP mass by measuring the energy spectrum of the recoiling 
nuclear targets. 86 Typically, of order 100 or more events are needed for such 
a determination to 10-20%. For higher WIMP masses, the recoil energy 
spectrum varies little, and WIMP mass extraction is much more difficult. 
Since the energy transfer from the WIMP to a nucleus is maximized when 
the two have the same mass, DD experiments with several target nuclei 
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ranging over a wide range of masses would facilitate the distinction between 
somewhat light and relatively heavy WIMPs, and so, potentially serve to 
establish the existence of multiple WIMP components in our halo. 

5.3. Indirect detection of neutralinos 

There are also a number of indirect WIMP search techniques that attempt 
to detect the decay products from WIMP annihilation at either the center 
of the sun, at the galactic center, or within the galactic halo. 

5.3.1. Neutrino telescopes 

Neutrino telescopes such as ANTARES or IceCube can search for high en- 
ergy neutrinos produced from WIMP- WIMP annihilation into SM particles 
in the core of the sun (or possibly the earth). The technique involves de- 
tection of multi-tens of GeV muons produced by interactions with polar 
ice (IceCube) or ocean water (ANTARES). The muons travel at speeds 
greater than the speed of light in the medium, thus leaving a tell-tale signal 
of Cerenkov light which is picked up by arrays of phototubes. The IceCube 
experiment, currently being deployed at the south pole, will monitor a cubic 
kilometer of ice in search of v M — > \i conversions. It should be fully deployed 
by 2011. The experiment is mainly sensitive to muons with > 50 GeV. 

In the case of neutralinos of SUSY, mixed higgsino dark matter 
(MHDM) has a large (spin-dependent) cross-section to scatter from hydro- 
gen nuclei via Z-exchange and so is readily captured. Thus, in the HB/FP 
region of mSUGRA, or in other SUSY models with MHDM, we expect ob- 
servable levels of signal exceeding 40 events/km 2 /yr with > 50 GeV. 
For the mSUGRA model, the IceCube signal region is shown beneath the 
magenta contour labelled fj, in Fig. 4. 87 These results were obtained using 
the Isajet-DarkSUSY interface. 71 Notice that DD signals are also observ- 
able in much the same region (below the contour labelled DD) where the 
neutralino is MHDM. 

5.3.2. Anti-matter from WIMP halo annihilations 

WIMP annihilation in the galactic halo offers a different possibility for 
indirect DM searches. Halo WIMPs annihilate equally to matter and anti- 
matter, so the rare presence of high energy anti-matter in cosmic ray events 
- positrons e+, anti-protons p, or even anti-deuterons D - offer possible 
signatures. Positrons produced in WIMP annihilations must originate rela- 
tively close by, or else they will find cosmic electrons to annihilate against, 
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or lose energy via bremsstrahlung. Anti-protons and anti-deuterons could 
originate further from us because, being heavier, they are deflected less and 
so lose much less energy. The expected signal rate depends on the WIMP 
annihilation rate into anti-matter, the model for the propogation of the 
anti-matter from its point of origin to the earth, and finally on the as- 
sumed profile of the dark matter in the galactic halo. Several possible halo 
density profiles are shown in Fig. 10. We see that while the local WIMP 
density is inferred to a factor of ~ 2-3 (we are at about 8 kpc from the 
Galactic center), the DM density at the galactic center is highly model- 
dependent close to the core. Since the ID signal should scale as the square 
of the WIMP density at the source, positron signals will be uncertain by a 
factor of a few with somewhat larger uncertainty for p and D signals that 
originate further away. Anti-particle propagation through the not so well 
known magnetic field leads to an additional uncertainty in the predictions. 
The recently launched Pamela space-based anti-matter telescope can look 
for e + or p events while the balloon-borne GAPS experiment will be de- 
signed to search for anti-deuterons. Anti-matter signals tend to be largest 
in the case of SUSY models with MHDM or when neutralinos annihilate 
through the A-resonance. 88 



10 n , 1 , 1 , 1 10" 




m (GeV) rn (GeV) 



Fig. 10. Various predictions for the DM halo in the Milky Way as a function of 
distance from the galactic center. The earth is located at r ~ 8 kpc. 
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5.3.3. Gamma rays from WIMP halo annihilations 

As mentioned in the Introduction, high energy gamma rays from WIMP 
annihilation offer some advantages over the signal from charged antiparti- 
cles. Gamma rays would point to the source, and would degrade much less 
in energy during their journey to us. This offers the possibility of the line 
signal from Z\Z\ — ► 77 processes that occur via box and triangle diagrams. 
While this reaction is loop-suppressed, it yields monoenergetic photons with 
Ej ~ mwiMP, and so can provide a measure of the WIMP mass. Another 
possibility is to look for continuum gamma rays from WIMP annihilation 
to hadrons where, for instance, the gamma is the result of it decays. Since 
the halo WIMPS are essentially at rest, we expect a diffuse spectrum of 
gamma rays, but with E 1 < towimp- Because gamma rays can traverse 
large distances, a good place to look at is the galactic center, where the 
WIMP density (see Fig. 10) is expected to be very high. Unfortunately, 
the density at the core is also very uncertain, making predictions for the 
gamma ray flux uncertain by as much as four orders of magnitude. Indeed, 
detection of WIMP halo signals may serve to provide information about 
the DM distribution in our galaxy. 

Anomalies have been reported in the cosmic gamma ray spectrum. In 
one example, the Egret experiment 89 sees an excess of gamma rays with 
E 1 > 1 GeV. Explanations for the Egret GeV anomaly range from Z1Z1 — > 
bb — > 7 with ~ 60 GeV, 90 to mis-calibration of the Egret calorimeter. 91 
The GLAST gamma ray observatory is scheduled for lift-off in 2008 and 
should help resolve this issue, as will the upcoming LHC searches. 92 



5.4. Gravitino dark matter 

In gravity-mediated SUSY breaking models, gravitinos typically have weak 
scale masses and, because they only have tiny gravitational couplings, are 
usually assumed to be irrelevant for particle physics phenomenology. Cos- 
mological considerations, however, lead to the gravitino problem, wherein 
overproduction of gravitinos, followed by their late decays into SM parti- 
cles, can disrupt the successful predictions of Big Bang nucleosynthesis. The 
gravitino problem can be overcome by choosing an appropriate range for 
m,Q and a low enough re-heat temperature for the universe after inflation 93 
as illustrated in Fig. 11, or by hypothesizing that the G is in fact the stable 
LSP, and thus constitutes the DM. 94 

Here, we consider the consequences of a gravitino LSP in SUGRA mod- 
els. If gravitinos are produced in the pre-inflation epoch, then their number 
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Gravitino Mass (GeV) 

Fig. 11. An illustration of constraints from Big Bang nucleosynthesis which re- 
quire Tfi to be below the various curves, for the HB/FP region of the mSUGRA 
model with mo = 2397 GeV, m 1 / 2 = 300 GeV, A$ — and tan/3 = 30, from 
Kohri et al. 93 to which we refer the reader for more details. 

density will be diluted away during inflation. After the universe inflates, it 
enters a re-heating period wherein all particles can be thermally produced. 
However, the couplings of the gravitino are so weak that though gravitinos 
can be produced by the particles that do partake of thermal equilibrium, 
gravitinos themselves never attain thermal equilibrium: indeed their den- 
sity is so low that gravitino annihilation processes can be neglected in the 
calculation of their relic density. The thermal production (TP) of gravitinos 
in the early universe has been calculated, and including EW contributions, 
is given by the approximate expression (valid for -C Mi 95 ): 

(^)(r&)"(io&v) <-) 

where Tr is the re-heat temperature. 

Gravitinos can also be produced by decay of the next-to-lightcst SUSY 
particle, the NLSP. In the case of a long-lived neutralino NLSP, the neu- 
tralinos will be produced as usual with a thermal relic abundance in the 
early universe. Later, they will each decay as Z\ — > 7G, ZG or KG. The 
total relic abundance is then 

n^h 2 = qTO + h 2 . (w) 
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The G from NLSP decay may constitute warm/hot dark matter depending 
in the Zi — G mass gap, while the thermally produced G will be CDM. 96 

The lifetime for neutralino decay to the photon and a gravitino is given 
by, 97 

48ttM 2 , r 2 
100 GevV 1 



r(Zi -> 7 G) A 2 - 

m| (1 - H) d (l + 3r 2 ) 



5.8 x 10 8 s 



m Si J A 2 (1 - r 2 ) 3 (l + 3r 2 ) 



where A = (t>£^ cosfl^ + v^ 1 ' sin9 w ) 1 , with Wg 1 ] being the wino and bino 
components of the Zi, & Mp is the reduced Planck mass, and r = wig/m^. 
Similar formulae (with different mixing angle and r-dependence) hold for 
decays to the gravitino plus aZor/i boson. We see that - except when 
the gravitino is very much lighter than the neutralino as may be the case 
in GMSB models with a low SUSY breaking scale - the NLSP decays well 
after Big Bang nucleosynthesis. Such decays would inject high energy gam- 
mas and/or hadrons into the cosmic soup post-nucleosynthesis, which could 
break up the nuclei, thus conflicting with the successful BBN predictions of 
Big Bang cosmology. For this reason, gravitino LSP scenarios usually favor 
a stau NLSP, since the BBN constraints in this case are much weaker. 

Finally, we remark here upon the interesting interplay of baryogenesis 
via leptogenesis with the nature of the LSP and NLSP. For successful ther- 
mal leptogenesis to take place, it is found that the reheat temperature of 
the universe must exceed ~ 10 10 GeV. 98 If this is so, then gravitinos would 
be produced thermally with a huge abundance, and then decay late, de- 
stroying BBN predictions. For this reason, some adherents of leptogenesis 
tend to favor scenarios with a gravitino LSP, but with a stau NLSP." 



5.5. Axion/axino dark matter 

If we adopt the MSSM as the effective theory below Mqut, and then seek 
to solve the strong CP problem via the Peccei-Quinn solution, 100 we must 
introduce not only an axion but also a spin-^ axino a into the theory. The 
axino mass is found to be in the range of keV-GeV, 101 but its coupling is 
suppressed by the Peccei-Quinn breaking scale f a , which is usually taken 
to be of order 10 9 — 10 12 GeV: thus, the axino interacts more weakly than 
a WIMP, but not as weakly as a gravitino. Both the axion and axino can 
be compelling choices for DM in the universe. 102 
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Like the gravitino, the axino will likely not be in thermal equilibrium 
in the early universe, but can still be produced thermally via particle scat- 
tering. The thermal production abundance is given by 102,103 

nTP , 2 cc6l /1.108\ /10 11 GeV\ 2 

/ rrig \ ( T R \ 
X V100 McV / V10 4 GeV / ' 

where f a is the PQ scale, N is a model-dependent color anomaly factor that 
enters only as f a /N , and g s is the strong coupling at the reheating scale. 

Also like the gravitino, the axino can be produced non-thermally by 
NLSP decays, where the NLSP abundance is given by the standard relic 
density calculation. Thus, 

n a h 2 = nT p h 2 + -^^n^ SP h 2 . (15) 

m N LSP 

In this case, the thermally produced axinos will be CDM for m 5 ~ 0.1 
McV, 102 while the axinos produced in NLSP decay will constitute hot /warm 
DM. 96 Since the PQ scale is considerably lower than the Planck scale, the 
lifetime for decays such as Z\ — > 7a are of order <~ 0.01 — 1 sec- well before 
BBN. Thus, the axino DM scenario is much less constrained than gravitino 
DM. 

Note also that if axinos are the CDM of the universe, then models with 
very large f^/i 2 ~ 100—1000 can be readily accommodated, since there is a 
huge reduction in relic density upon Z\ decay to the axino. This possibility 
occurs in models with multi-TeV scalars (and hence a multi-TeV gravitino) 
and a bino-like Z\ . In this case- with very large rriQ- there is no gravitino 
problem as long as the re-heat temperature Tr ~ 10 6 — 10 s GeV. This range 
of Tr is also what is needed to obtain successful non-thermal leptogenesis 
(involving heavy neutrino N production via infiaton decay) 104 along with 
the correct abundance of axino dark matter. 108 

5.5.1. Yukawa-unified SUSY with mixed axion/axino dark matter 

The gauge group 5*0(10) is very highly motivated in that it unifies all 
matter particles of each generation into a spinorial 16 dimensional rep- 
resentation. Included in the 16 is a state containing a SM gauge singlet 
right-hand neutrino. SUSY 5O(10) theories may also allow for third gen- 
eration Yukawa coupling unification. In Ref., 105,106 it was found via scans 
over SUSY parameter space that in fact it is possible to find models with 
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Yukawa coupling unification, but only for certain choices of GUT scale SSB 
boundary conditions. Yukawa unified solutions can be found if 

• mi6 ~ 3 - 15 TeV, 

• mi/2 is very small, 

• tan/3 ~ 50, 

• Ao ~ — 2mi6 with mio ~ 1.2mi6, 

• mjj d > m ffu at M GUT . 

In this case, models with an inverted scalar mass hierarchy are found. The 
spectra is characterized by 1. first/second generation scalars in the 3-15 
TeV range, 2. third generation scalars, tua and ji in the few TeV range 
and 3. gluinos around 350-500 GeV, with charginos <~ 100 — 160 GeV and 
Zi ~ 50 - 90 GeV. 

A problem emerges in that the calculated ncutralino relic density is 
about 10 2 — 10 5 times its measured value. A solution has been invoked that 
the axino and axion are instead the DM particles. In this case, neutralinos 
would still be produced at large rates in the early universe, but each neu- 
tralino would decay after a fraction of a second (slightly before the onset 
of BBN) to axino a plus photon. Then the relic axino density would be 
(mj/ro^JO^/i 2 , and the mass ratio out in front reduces the relic density 
by a factor of 10~ 2 — 10 -5 !. The neutralino is still long lived enough that 
it will give rise to missing energy at the LHC. In fact, since gluinos are 
so light, we would expect a gluino pair cross section of order 10 5 fb, along 
with decays g — > 66Z 2 , bbZ\ and tbW\ + c.c.. 107 These new physics reac- 
tions should be easily seen via isolated multi-muon plus jet production in 
the early stages of LHC running. 

The Yukawa-unified SUSY scenario is also very appealing cosmologi- 
cally. We would expect the gravitino mass to be of order m\§: 3-15 TeV. 
If it is heavier than about 5 TeV, then the gravitino decays with lifetime 
around 1 sec: right at the onset of BBN! This eliminates the BBN gravitino 
problem as long as re-heat temperature Tr ~ 10 9 GeV. Several scenarios of 
mixed axion/axino dark matter in Yukawa unified SUSY have been exam- 
ined in Refs. 108 ' 109 It is found that one may accommodate either dominant 
axion or axino DM, but mi6 at the high end. Values of m\e ~ 10 TeV are 
preferred, as are high values for the PQ breaking scale: f a ~ 10 12 GeV. The 
values of Tr allowed can accommodate several baryogenesis mechanisms, 
for instance, non-thermal leptogencsis. 
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6. Conclusions 

The union of particle physics, astrophysics and cosmology has reached an 
unprecedented stage. Today we are certain that the bulk of the matter in 
the universe is non-luminous, not made of any of the known particles, but 
instead made of one or more new physics particles that do not appear in the 
SM. And though we know just how much of this unknown dark matter there 
is, we have no idea what it is. Today, many theoretical speculations seek 
to answer one of the most pressing particle physics puzzles, "What is the 
origin of EWSB and how can we embed this into a unified theory of particle 
interactions?" The answer may automatically also point to a resolution 
of this 75 year old puzzle as to what the dominant matter component of 
our universe might be. Particle physicists have made many provocative 
suggestions for the origin of DM, including supersymmetry and extra spatial 
dimensions, ideas that will completely change the scientific paradigm if they 
prove to be right. 

The exciting thing is that many of these speculations will be directly 
tested by a variety of particle physics experiments along with astrophysical 
and cosmological searches. The Large Hadron Collider, scheduled to com- 
mence operation in 2009, will directly study particle interactions at a scale 
of 1 TeV where new matter states are anticated to exist for sound theoreti- 
cal reasons. These new states may well be connected to the DM sector, and 
so in this way the LHC can make crucial contributions to not only particle 
physics, but also to cosmology. If indeed the LHC can make DM particles 
or their associated new physics states, then a large rate for signal events 
with jets, leptons and $t is expected. 

Any discovery at LHC of new particles at the TeV scale will make a 
compelling case for the construction of a lepton collider to study the prop- 
erties of these particles in detail and to elucidate the underlying physics. 
Complementary to the LHC, there are a variety of searches for signals from 
relic dark matter particles either locally or dispersed throughout the galac- 
tic halo. The truly unprecedented thing about this program is that if our 
ideas connecting DM and the question of EWSB are correct, measurements 
of the properties of new particles produced at the LHC (possibly comple- 
mented by measurements at an electron-positron linear collider) may allow 
us to independently infer just how much DM there is in the universe, and 
quantitatively predict what other searches for DM should find. 
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